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^^ • Abstract 



Scattering of a plane electromagnetic wave by an anisotropic impedance concave wedge at skew 
incidence is analyzed. In the case of a right-angled wedge, the problem reduces to two symmetric 

O ; order-2 vector Riemann-Hilbert problems (RHPs) which are solved exactly. The problem of 

matrix factorization leads to a scalar RHP on a genus-3 Riemann surface. Its closed-form solution 
is derived by the Weierstrass integrals, and the associated Jacobi inversion problem is solved in 

C^ I terms of elliptic integrals. It is shown that the solvability of the physical problem is governed by 

the location of the four zeros of the product of two quadratic polynomials associated with the 
two vector RHPs. The coefficients of these polynomials are expressed through the wave number, 
the impedance parameters and the angle of incidence. In the case when the solution is unique, 

^ I the reflected, surface and diffracted waves are recovered. 

00 ' 

0\ . 1 Introduction 

en 

^J . In diffraction theory the study of scattering of a plane obliquely incident electromagnetic wave from 

an anisotropic impedance wedge is one of the key canonical problems. In the case of normal incidence 
when the tensor impedance has zero diagonal and nonzero off-diagonal entries, the problem was 
solved in ([1]). For oblique incidence, a closed-form solution is known for some special cases including 
the vector case for a half-plane ([2]), ([3]), ([4]), the diagonal case for a wedge when the impedance 
parameters meet certain conditions, and the problem reduces to decoupled Maliuzhinets equations 
C^ ' (Ej) , and the triangular case for a wedge when the impedance tensor allows for the reduction of the 

problem to two consequently solved scalar Maliuzhinets equations ([6]). 

Approximate solutions by the perturbation technique are available for the cases when the wedge 
is almost a half-plane ([7j) or when the incidence is almost normal ([8]). Different approximate nu- 
merical solutions are available for the general case. They include those obtained by the method of 
integral equations ([9]), (1101) . the method of approximate matrix factorization in conjunction with 
the Fredholm integral equation theory (jlip . and the probabilistic random walk method (|12|) . A 
method based on splitting the spectra into two functions, the solution of the Maliuzhinets equation 
and a function that is defined by a series whose coefficients solve a certain infinite system of linear 
algebraic equations, was recently proposed in (I13J) . By this method, the diffraction field in the case 
of scattering by the wedge of angle 37r/2 was recovered in (J14p . A method of reduction of wedge 
diffraction problems to functional equations was presented in (I15p . In particular, it was pointed 
out that the electromagnetic problem for a right-angled convex wedge can be reduced to an order-6 
vector RHP. 



An exact solution for the general case of the tensor impedance and obliquely incident electromag- 
netic wave even for a right-angled wedge is still unavailable in the literature. Our goal in this paper 
is to derive a closed-form solution of the problem on an impenetrable right-angled concave wedge 
(of angle 37r/2) in the case of oblique incidence when on the faces S^ of the wedge the boundary 
conditions are 

±£, = rj%znp, T£p = vtzn,, {p,e)GS^. (i.i) 

Here, S^ = {9 = 0,0 < p < oo}, S^ = {9 = -ir/2,0 < p < oo}, Sp^l-Lp and £z,T~{-z are the p- and 
z-components of the electromagnetic field, Z is the intrinsic impedance of the medium, and t]^ and 
r/^ are the impedance parameters. 

In Section 2, we constitute the boundary value problem for the exterior of the concave wedge 
of an arbitrary angle, We = {0 < p < oo,0 < 9 < a} (0 < a < vr), with the general impedance 
boundary conditions 

±£, = r^%Z'Hp + r^%Zn., T£p = ritpZTip + r^tZTi,, ip,e)GS^. (1.2) 

Next, by applying the Laplace transform (jlSp we reduce the problem to two sets of six functional 
equations. We focus on the case a = 7r/2, and in section 3, we reduce the functional equations 
to two vector RHPs with a 2 x 2 matrix coefficient. The solution to the second RHP is derived 
from the solution to the first RHP by a certain transformation of the impedance parameters and 
the angle of incidence /? G (0, vr). The solution of the RHPl and RHP2, vectors $^(r/) and ^^(77), 
respectively, have to satisfy the symmetry conditions <^"'"(r/) = <1>~(— r/) and ^"'"(ry) = ^^(—rj). It 
is shown that the matrix coefficient has the structure G{rj) = hi(rf)Qi{r]) + b2{r])Q2{r]), where bi(r]) 
and b2{rj) are Holder functions on the contour of the problem L = {—00, +00}, and Qi{rj) and Q2{fl) 
are polynomial 2x2 matrices. This case is algebraic, and the symmetric vector RHPs can be solved 
exactly ([16]), m\\ . 

In section 4, for simplicity, we assume that the impedance parameters rjpz and r]zp vanish, that 
is the boundary conditions (|1.2p become (jl.ip . For this case, we reduce the problem of matrix 
factorization to a scalar RHP on a genus-3 Riemann surface. The coefficient of the RHP on the 
first and second sheets coincides with the first and second eigenvalues of the matrix coefficient, and 
depending on the impedance parameters and the angle of incidence, the index of the quotient of 
the eigenvalues may be 0, 1 and —1. Because of the symmetry of the problem we manage to solve 
the associated Jacobi inversion problem in terms of elliptic integrals, and the solution to the matrix 
factorization problem is constructed in closed form by quadratures. 

In section 5, we derive a closed-form solution to the vector RHPs. It turns out that the number of 
free constants in the general solution is governed by the location of zeros of a certain polynomial whose 
coefficients depend on the impedance parameters and the angle of incidence. The free parameters 
have to satisfy the symmetry conditions, additional mathematical and physical conditions. The 
results are summarized in Theorem 5.1 which implies that the solution is unique if and only if all 
the four zeros of the characteristic polynomial 

p(r?) = [(7?2 - kl) cos2 /3 - (7? - 7r)(r/ - ^^)][{r^^ - kl) cos^ /3 - (r? - 7i+)(r/ - 74+)] (1-3) 

lie in the upper half-plane. Here, 7^^ = k{r] p)~'^ siv? 13 , 74 = k'q'f^ siv? fi. In the case of normal 
incidence, /3 = '7r/2, a simple closed-form solution is obtained, and the characteristic polynomial p{rj) 
reduces to p{r]) = {rj — Ji){ri — Ji){ri — ^^){r] — jf). The solution is unique if Im7- > (j = 1, 4). 
This is consistent with the uniqueness theorem proved in (jlSp for an arbitrary wedge in the scalar 
case of normal incidence. 

In section 6, we derive the spectra Se{s) and Sh{s) of the problem and express the Sommerfeld 
integrals through the solution to the two vector RHPs. On continuing analytically the functions 
Se{s) and Sh{s) to the right and to the left, applying the steepest descent method and the Cauchy 
theorem, we find the reflected, surface and diffracted waves. 



2 Formulation: diffraction by a concave wedge of an arbitrary angle 

Consider diffraction of an electromagnetic wave 

(^gtnc^ j^inc^ ^ l^^i^ jji^^^ikz cos P-iut^ (2.1) 

where 



[eI, Hi) = [Eo, Ho) e~ 



ikpcos{9—9o) sin (9 



(2.2) 



by a wedge W = {0 < p < oo,a < 9 < 2tt, \z\ < oo}, a £ (0, vr), cliaracterized by tlie impedance 
boundary conditions 



±£z = ri^pZUp + ry^^ZUz on 5=^, 
=\^£p = 7]fpZnp + r]f^Znz onS^. 



(2.3) 



Here, k (Im A; > 0) is the wave number, uj is the angular velocity, Z = v/W^ ^^^ intrinsic impedance 
of the medium, po is the magnetic permeability, Eq is the electric permittivity, rjpp, rjpz, rjzp and rj^z are 
the impedance parameters. The symbols S^ stand for the boundaries of the wedge, S~ : 6 = 0, p > 0, 
and 5*+ : 9 = a, p > 0. Because of the representation ()2.ip for the incident waves it is natural to 
split the electric and magnetic fields as 



£, n 



E, H]e 



^ikz cos /B—iLut 



(2.4) 



By employing the Maxwell equations we eliminate the radial components of the electric and magnetic 
fields 



E, 



ZH, 



k sin"' /3 

i 
ksin j3 



d Z d 

cos/3—Ez + -—Hz 
op p oO 

Zcosl3—Hz--—Ez 
op p OO 



(2.5) 



The resulting boundary conditions, after rearrangement, are formulated in terms of the z-components, 
(pi = Ez and (/>2 = ZHz, as 



i dcpi 



,502 



I U(pi . o(p2 _^ ± , , ± i, n c± 

p 06 op 



where 



d0i id4>2 ± ± ± 

zcos/;^^ \ — - +73 01 ±74 02 = on 6 , 

op p off 

7i = ^, 72 = 



'Tw^' 



r^P 



± r?±A;sin p ^ ks\v? p ^ ^ ^ ^ 

l3 = H ' 74 = ^± iVppVzz - VpzVzp) 



^ 



^P 



(2.6) 



(2.7) 



In the exterior of the wedge. We = {O<p<oo,O<0< a}, the functions 0i and 02 satisfy 
the Helmholtz equation. To convert the boundary value problem for the Helmholtz equation to a 
system of functional equations, we use the scheme (I15p . On making the affine transformation of the 
Cartesian coordinates u = x — ycoia and v = ycsca we obtain that, in the new coordinates, the 
functions 0i and 02 are solutions of the differential equation in a quarter-plane 



7-^ + TT^ - 2 COS a——- + fco sm a 
ow^ ov-^ ouov I 



0, < u,v < OO, j = 1, 2, 



subject to the four boundary conditions 
1 d d 



cot a— j (pi+i cos /3— 7i'</'i + 72" 'A2 = 0, {u, v) G S~^ , 

icos/3— i -i {- cot a--) (t>2 + 73 (pi +lt<p2 = 0, {u,v) € 5+, 

ov V sm a ou ov J 

cota^- I (/>! +icos/3— 1- 7~(/);^ + ^"(^2 = 0, (n, w)g5~, 



sma 



dv du I du 



icosf3--^ + i[- cota-- ) 02 + 73"0i -7^02 = 0, {u,v) <G S' , (2.9) 

ou \sm.adv ouj 

where ko = ksm(3, 5+ = {{u,v) G M^ : u = 0,u > 0}, 5" = {{u,v) G R^ : f = 0,u > 0}. By means 

of the Laplace transform 

00 

4)j{u,r])= [ e''^''cl)j{u,v)dv (2.10) 



the differential equation (|2.8p may be put into the form 

\— + 2ir]cosa— + kls\v?a-rj^\(t)j{u,ri)=fj{u), j = 1,2, (2.11) 

where 

fj{u) = ^(u,0) - ir](l)j{u, 0) - 2 cos a^{u,0). (2.12) 

The two roots of the characteristic equation of the differential operator in (I2.1ip are —ir] cos a it 



r/^ — /cg sin a. Fix a branch C, = Jrf' — k^ of the two- valued function C,"^ = ff — A^q by the condition 

C(0) = —iko. The branch is a single-valued analytic function in the T/-plane cut along the line joining 

the branch points ib/co and passing through the infinite point. As —00 < r/ < -|-oo, the branch chosen 

possesses the property Re^ > 0, and the general solution to equation (12. lip bounded as n — )• 00 has 

the form 

cPjiu,v) = ^i(r?)e-«" - - / e-l"-"il«/,(ni)dni, j = 1,2, (2.13) 

2t, Jo 

where ^ = — ir/ cos a -|- C sin a. To derive a functional equation for cpj, we apply the Laplace transform 
with respect to u using the function ^ as its parameter 

^Ai^^v)= e-^''(l)Au,v)du. (2.14) 

JO 
On referring to (j2.12p . we find that 

/■OO fj 

/ fj{u)e~^''du = —(t)Aii,Q)-{ir] + 2icosa)(t)j{ii,Q) + 2cosa(l)°, (2.15) 

Jo dv ■' 

where (j)°. = 0j(0, 0). On the other hand, from (I2.13P 

1 poo 



UlfJj 

du 



2Uo 

1 f°° 
{0,rj) = -CA,irj)--J^ e'^^f,{u)du. (2.16) 



By excluding Aj{i]) and the Laplace transforms of the functions fj{u) and employing the second 
relation in (j2.16p we get the following two equations: 

^{0,7]) + C4>j{0,rj) + ^{iC,0) - {ir] + 2Ccosa)^j{i^,0) + 2cosa(l)° = j = 1,2. (2.17) 

These two equations should be complemented by the Laplace-transformed boundary conditions. We 
apply the Laplace transform (I2.10p to the boundary conditions (12. 9j) on the boundary S~^ and the 
transform (j2.14p to the conditions on the side S'. We have 

icsca— — (0,77) — (7^ + r/cot a)(^i(0, 77) + (7^ + ry cos /3) 02(0,??) + icot a0j — icos/302 = 0, 

— icsca— — (0,r/) + (73" + t/cos/3)i^i(0, 77) + (74^ + 77 cot 0)^2(0,?/) — i cos I3(pl — icot a(p2 = 0, 

— icsca— — (i^,0) + (7^ + i^cot a)</)i(z^,0) + (7^ + i^cos f3)(j)2{iS.,0) — icot acpl — icos I3(j)2 = 0, 

zcsccK— — (iC)0) + (73^ -|-i^cos/3)(/>i(i^,0) — (7^ -|-i^cot a)(/>2(iC,0) — z cos /3(/)° + i cot acf)^ = 0. (2.18) 

Notice that by applying the Laplace transform (I2.10p with respect to u and then the Laplace trans- 
form (j2.14p with respect to v we obtain another set of six equations which coincide with (|2.17p and 
(I2.18P if we make the following transformation of the parameters and the functions: 

li^^li^ it'^-li^ 73" ^ -73"; 74" ^71' /3^7r-/3, 
,^,(0,r?) o ^,(r?,0), ^(0,??) ^ ^(^>0), 

(^j(ie,0)f^(^,(0,ie), ^(^^'0)^^(0'^^)' i = l'2. (2.19) 

3 Vector RHPs in the case of a right-angled concave wedge 

Consider the case a = 7r/2 (Fig. 1) and reduce the system of six functional equations (I2.17|) . (I2.18P 
and the one obtained from the system (|2.17p . (|2.18p by the transformation (j2.19p to two symmetric 
RHPs with a 2 X 2 matrix coefficient. If the domain of interest is a quarter-plane, then the standard 



Cartesian coordinates x = u, y = v can be employed, and ^ = C^ C = yv'^ ~ ^o- ^^ ^P^^* ^^^ incident 
waves El and HI into two waves 

i^El+, Hi+) = [Eo, i^o) e-*"'^°^(^-^«)a;(0;eo, ^2), 

(^r, HI-) = [Eo, /7o)e-*'=»^^°^(^-^°)u;(0;O,eo), (3.1) 

where 

Consider first the problem associated with the incident wave {El'^,Hl'^). From the boundary condi- 
tions (|2.18p we express the derivatives d(j)j{Q,rj)/dx and d(j)j{i(,0)/dy through the functions (l)j{0,r]) 
and (pj{i(,0) {j = 1,2). Then, on replacing r] by —r] in (|2.17p we obtain two more equations 

^(0,-77) + C<^,(0,-?7) + ^(iC,0)+i7y,^,(iC,0) = 0, i = l,2. (3.3) 



incident wave 



point of observation 




Figure 1: Geometry of the problem on a right-angled concave wedge. 



These two extra equations in conjunction with (12.171) allow us to express the functions (pj{iC,0) 
through the functions (pj{0,±r]) 

M<, 0) = ^ [</'i(0, W - M^, -r])] H ^ 02(0, rj) + ^ '^'2(0, -r?), 



2t] 



27? 



2?? 



M^CO) = -^^±^^,(0,,) + ^^-r^^ MO, -V) - ^^[<^2(0,,) - 0.(0,-,)]. (3.4) 



2r/ 



2r? 



2r/ 



On substituting the expressions for the functions (f)j{iC,,Q) and the derivatives d(j)j{0,r])/du and 



d(j)j{iC,0)/dv in equations (j2.17p we ultimately obtain 

^(77)$+ (r/) + B{rj)<^+{-r]) =0, r] e L = (-00, +00), 
where 



(3.5) 



$+(r?) 



0i(O,r/) 



^(r/) 



aii(r/) ai2(??)\ ^. ^^ ^ /&ii(??) &12(??) \ 



02(0, r/)^' ^^^''^ \^ 021 (r/) 022(^7)/ ' \b21iv) b22i'n) J 

an = 5i+(r/)[l - 5i-(r/)] - (52-(??)(^3+(r/), ai2 = -52+(r?)[l - (5i-(r?)] - 62-{v)h+{v), 

a2i = S3+{r])[l - (54- (r?)] + d3-{r])6i+{r]), 022 = '54+(r/)[l - 54-(r/)] - 63-{v)S2+{v), 

bu = 6i-{r])6i+{-r]) + 62-{r])6-i+{-r]), 612 = -(5i_(r/)(52+(-??) + 62-{v)S4+{-r]), 

b2i = -(53_(r/)(5i+(-r?) + 64-{v)S3+{-r]), 622 = (53_(?7)(52+(-r?) + 54_(r/)(54+(-r?), 



^i+(^)=7l+<, 5i-(??) 



^i+ iv) =lt +^ cos /3, (5j_ (r?) - ^ 



2i +^ 
7~ + iC, cos /? 



, j = l,4. 



2?? 



, j = 2,3. 



(3.6) 



It can be directly verified that A{ri) = B{—rj) for all rj. Denote $ (77) = $+(—77) and let r/o = ko sin^o- 
Then the vectors <l>^(?7) constitute the solution to the following symmetric RHP. 



RHP 1. Find two vectors <^^(r/) analytic everywhere in the half-planes C^ : ztlmry > except 
at the simple poles rj = ±770, 

res ^tiri) = iEo, res ^tiv) = iZHo, (3.7) 

and continuous up the the real axis L. At the contour L, their limit values satisfy the boundary 
condition 

<^+{7]) = G{rj)<^-{r]), r]GL, (3.8) 

where G{rj) = —[A{rj)]^^B{rj). At infinity, both vectors vanish, and everywhere in the plane they 
meet the symmetry condition ^~^(r]) = ^^{—t]). 

The same derivations can be applied to the case of the incident wave {El" , HI" ) . On making the 
transformations (j2.19p in the system ()2.17p . ()2.18p we arrive at the second RHP. 

RHP 2. Find two vectors ^ (77) analytic everywhere in the half-planes C except at the simple 
poles T] = ±fjo, 770 = fco cos 6*0, 

res *j^(7?) = iEo, res ^^(r/) = iZHo, (3.9) 

V=VO V=VO 

and continuous up the real axis L. At the contour L, their limit values satisfy the boundary condition 

^+{7])=G{7])'^-{7]), 7]£L, (3.10) 

where ^^(77) = 0^(77,0), "^^{7]) = ^H— 77), j = 1,2, and G{i]) is the matrix G{r]) whose parameters 
7- (j = l,...,4j and /3 are transformed by i2.19\) . At infinity, both vectors, "^^{rj) and ^"(77), 
vanisii. 

It is sufficient to solve the RHP 1. The solution to the RHP 2 can be derived by the transformation 
(12.19p . Notice that the matrix coefficient of the RHP admits the representation 

Giv) = ^[G'{v)+CGHv)], (3.11) 

where ^(77) is a scalar, G^ and G^ are 2x2 polynomial matrices, and det G{r]) 7^ 0, r/ G L. It is known 
(Antipov, 2010) that any nonsingular matrix of such a structure can be factorized by reducing the 
vector RHP to a scalar RHP on a Riemann surface of the associated algebraic function. To simplify 
our derivations, it will be helpful to assume that either 72 =73 =0 and the angle /3 is arbitrary, or 
that all the four impedance parameters 7- are arbitrary and /3 = 7r/2. In what follows we consider 
the former case, namely, 72 =73 = 0, /? G (0,vr). 

Before beginning factorization of the matrix G{r]), we notice that the matrix ^(77) is continuous 
everywhere in the contour L, and as r/ — t- ±00, G(r/) ~ diag{— 1, —1}. Compute first the index of its 
determinant 

Ko = T^[argdetG(r/)]-^, (3.12) 



which has the form 



and 



27r' 
•5o('f) = ('f - ko) ™' P-(v--ti)(l- 7i\ (3.14) 



We distinguish the following three cases: 

(i) the two zeros of the quadratic polynomial 6q{i]) lie in the lower half-plane, (5o(^) = ~(^ + 
n){ri + ra) sin^ /?, Imr,- > 0, j = 1, 2, 



(ii) the zeros lie in the opposite half-planes, 5o(^) = ~{v ~ Ti)iv + ^2) sin^ /?, Im'^j > 0, j = 1, 2, 
and 

(iii) both zeros lie in the upper half-plane, 6o{ri) = —{rj — Ti){r) — T2) sin^ /3, Imrj > 0, j = 1, 2. 

By the argument principle, in case (i), kq = 2 (in the upper half-plane, the function So{—ri) has 
two zeros, while the function 6o{r]) has no zeros in C"^). In case (ii), kq = 0, and in the last case, 
Ko = -2. 

Our next step is to rearrange the representation (j3.1ip as 



G(7?) 



c 



:GH^) 



1 + 



1 



ri{ri) r2(r/) 



^o(r/)(5i(r?) 
where / is the unit 2x2 matrix, rj{rj) (j = 1, . . . , 4) are polynomials, 

Siiv) = (71" + iOilt + iC) + V^ cos^ /3, 
the entries gltii]) of the polynomial matrix G^{r]) are 

9lj(^) = -^(Ti" + 74")[^o cos^ /3 - r/2 sin^ /3 + (-l)^'r/(7f - 74") + 71-74" 



(3.15) 



(3.16) 



1,2, 



Suiv) = -2irycos/3[7i 74 - (74")^ - fc^sin^/?], 

5^1 (r?) = 2ir?cos/3[7f 74- - (7+)^ - k^ sin^ /3], (3.17) 

and di(??) = det G^{r]) is a degree-4 even polynomial 

di(ry) = Qo + airy^ + 02^?^, (3.18) 

whose coefficients are 

ao = -(7i+ + 74+)^(7r74" + ^0 cos^ /3f, 02 = -(71+ + 74+)^ sin^ P, 

«i = (7i+ + 74+)'[(7r)' + (74-)' - 2feo' cos^ /3] + {-4(7r74-)' " (27i+74+ " k'of 
+ 271-74" [(7i+ - 74+)' + 2A:o'] + 2A:2[(7i+)' - 271-74" + (74+)' + ^0] cos 2/3 - A;^ cos^ 2/3} cos^ /3. (3.19) 
It is an easy matter now to represent the matrix G(r/) in the form 



a(„) = '-i^aHmn}. 

where 5^ is a constant to be determined, 

bo{r]) + Co{r])l{r]) co{v)m{r]) 



(3.20) 



r(r/) 

6 



Co 



Co{ri)n{ri) bo{r]) - Co(r?)/(r/) y ' 

= -^, A = b^-c^f, f = f+mn, 

6 = ^-(^l + — ^ ^^i7??cos^ 
(5*(5i V C^i/ ' S^:6idi 

- - I + + ; 2 ■ 2 o 

7 = 7i 74 + 7i 74 - «o sm /3, 
and I, m, n and r are polynomials, 

/(??) = r/cos/3[r?2(7j^ - 7^) - 2r/(7f - 7^)(7i+ + 74^) 



(3.21) 



+(7i+ - 7|)(27r74" - 3A;g) - (7i+ - 74+)(^' " ^o) cos 2/3], 
Hv) = -^v^ + vHli - ll) + Mil - l4)[(.ltf - ^o] 

+ [(74+)' - A:o'](27r74- + ^o) + |^'[47r74" " 4(74+)' + k^] - \v\v' " ^o) cos4/3 

+ {V^ - r]\l^ - ll) + Alill - {it? - kl] + kl[{-itf - kl]} cos 2/3, (3.22) 

the polynomial n{ri) coincides with —m{rj) if 7j~ and ^f are replaced by 7^ and 74", respectively, 
and 

r{v) = ^{it + 74'') [-^' + 2(71+74+ + ^o) + rl" cos 2/3] 
x[3r/4 - Srfi^^f - Srf^^^l - 8r]\j^f + 8(7r)'(74")' + ^v^k^ 
+871" 74" ^0 + 3A:^ - 4(j?2 _ A;2)(r,2 ^ 2^-7" + k^) cos 2/3 + (t?^ _ ^2)2 cos 4/3]. (3.23) 

In ()3.2ip we divided the coefficients of the matrix T[r]) by ^/A[r]) (•\/A(7/) is a fixed branch of the 
function A^''^{r])) and in (j3.20p we multiplied the matrix r{7]) by the same factor. This transformation 
implies detr(r7) = 1 for all 77 € L and simplifies the solution. Notice that from ()3.20p . 

On comparing (|3.24p and (j3.13p . we obtain that A{ri) is a rational even function 

A,„ . «^. (3.5) 

By the argument principle, its index is equal to 0. 

4 Scalar RHP on a Riemann surface 

4.1 Statement of the problem 

Since G^{r]) is a polynomial matrix, the problem of factorization of the matrix G{rj) reduces to the 
one for the matrix T{r]). Its factorization can be expressed through the solution of the associated 
scalar RHP on the genus-3 Riemann surface 91 of the algebraic function ■w'^ = f{rj). The function 
/(??) is an even degree-8 polynomial which has the form 

/(r/) = /lo + hiri" + /i2r?^ + h^rf + h^rj^ (4.1) 

where 

ho = —4:UiU2U^, hi = —4 sin /3, 

hi = 4(7f - 74')'^^iW2 + (7+ - 74+)'(27r74" " 3A:g + kl cos 2(3f cos^ /3 - 2^3 

x{n2[(7r74~ - (7i+)')(cos 2/3 + 3) + 2kl sin^ /3] + ui[(7r74" " (74+)')(cos 2/3 + 3) + 2A;2 sin^ /3]}, 

/i2 = 'iulu^ + 8(7^747^4 + 7]^74^n4 ) cos^ /3 + 87^7]^ 7^74" cos^ /3(cos 2/3 - 3) 

-[(7r74~)' + (71+74+)'] (cos 2/3 + 3)^ + 2A;2 sin^ /3{-n+(3cos 2/3 + 1) 

+2 cos^ /3 [-271+74!" (cos 2;3 - 3) + 7i"7^(cos 2/3 + 3)] - 4nJ - 2A:g sin^ /3(1 + cos'^ /?)}, 

/i3 = 4sin^ /3[uJ - n+ + 2(7^74" - 7+7+) cos^ /3 + 2A:g sin^ /3], 



m = (7+)2 - /eg, U2 = {jtf-kl n3 = 71-74- + A:2cos2/3, uf = (iff + {jf? ■ (4.2) 



Let zboj (j = 1,2,3,4) be the eight branch pomt of the function \/f(j]) (they are determined 
expUcitly by radicals) such that < Imai < Ima2 < Ima^ < Inia4. Cut the plane along the 
segments l^ and I2 with the starting points ztoi and ±03 and the terminal points ±02 and ±04 , 
respectively. Fix the single branch of the function y/fjj]) as follows: 



/(ry) = 2i sin^ /3 J] ./ptpJe^^'^J^/' , (4.3) 

i=i 

where pf = |r/=Faj|, iff = arg{r]^aj), (^+ € (-vr,0), ipj G (0,7r), r] e L, j = 1,2,3,4. 

The genus-3 hyperelliptic surface IH is formed by two copies, Ci and C2, of the extended ry-plane 
C U {00} cut along the segments l^ and I2 ■ The two sheets are glued according to the rule 



^_1 ^/m, ^eci. 



Let \{rj, w) be a new function defined on the surface y\ as 



A(^>^) = I ^S'^l' ''''?' (4-5) 

^'' ^ I A2(??), r/GC2, ^ ^ 



where Ai(7?) = bo{r]) + co{r])\/ f{rj) and Ai(r7) = 6o(^) — co(r?)-\//(??) are the eigenvalues of the matrix 
r(7/). We find it convenient to introduce two matrices, Y and Q, as 



^(^>f^) = 2 



/+-Q(r/) 



^^'^Uw't) ■ <«' 



Then a solution to the matrix factorization problem 

r(r?)=X+(77)[X-(7?)]-\ r?GL, (4.7) 

may be expressed through a nontrivial solution to the following scalar RHP on the surface D\. 

Find a function x(??, w) piece-wise analytic on the surface *H\ £, £ = (L C Ci) U (L C C2), except 
for at most a finite number of poles, continuous up to the contour £, bounded at both infinite points 
of the surface ?l and satisfying the boundary condition 

X^{r],w) = \{ri,w)x'{r],w), {rj,w) e £.. (4.8) 

In terms of the solution of the RHP ()4.8p . the matrix X{ri) and its inverse are given by 

Xir]) = x(r?, w)Y{7], w) + xiv, -w)Y{r], -w), 7] G C^. 

[X(,r' = ^ + ^^^. .€C±. (4.9) 
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4.2 Solution with an essential singularity at infinity 

The chief difficulty in the procedure of matrix factorization based on the use of the solution of 
the scalar RHP on a Riemann surface (j4.8p arises in the necessity of constructing a meromorphic 
solution with explicitly determined poles and zeros. In order to derive such a solution, we consider 
the function 

Xo{r],w) =exp{'>po{7],w)}, (r/, u;) G <K, (4.10) 

where ipQ{r]^w) is the Weierstrass integral 

V^o(r/,u^) = ^ /logA(t,e)^-^, (4.11) 

Ml Js, IE, t — rj 

and ^ = w{t), t £ £,. This integral over the contour C on the surface IH can be transformed into two 
integrals over the contour L in the //-plane as 

Mv,w) = ±-,j[\og\,{t) + \og\2{t)]^^ + £.J^[iogX,{t) - logA2(t)]^=^-^. (4.12) 

To simplify this representation, we study the eigenvalues Xj of the matrix r(?7). Referring to (j3.2ip . 
we obtain that 6(r/) = b{—ri), c{r]) = —c[—r]) and also 



As ?7 — )• oo. 



A,MA,h)=detr(„) = l, ^ = M±£Mi^ = MZ2) , (4.13) 

A2(?/) K??)-c(r/)v/M Ai(-?7) 



7 ^ COS [) 

^* ill +7I) Till +1I) sm" /? 



On putting 6^, = —iijf + 7^) ^, we have 6(77) = 1 + 0{ri ^), rj — ;■ 00. At zero and at infinity, the 
function Xi{r]) / X2{r]) has the following properties: 

It has also been established that Xi{rj) and A2(^) are bounded and do not vanish on the contour L. 
On fixing the branches of the functions log Ai(7y) and X2{r]) as logAi(oo) = logA2(oo) = we obtain 
that logA2(??) = — logAi(7/) for all 7/ G L and therefore 

MvM = —[ r^^'^' ■ (4.16) 

where 

£(r?)=log/(r?), /(r?) = ^, (4.17) 

and the branch of the logarithmic function is fixed by the condition e(oo) = 0. To establish whether 
the function e(r/) is continuous in the contour L, we study next the increment of the argument of 
the function /(r/) when r/ traverses the contour L in the positive direction. Introduce the index of 
the function l{rj) 

ind/(77) = ^[arg/(77)]~^. (4.18) 
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Figure 2: Parametrically defined function /(?y), < rj < +00, in case (i, section 3) wlien ko = 1 + 0. li, 
/3 = 7r/4. Case a: kq = {'yf = 1 — i, 7^ = 1 + 2i, 7{" = 1 — 2i, 7^ = 1 — 3i). Case b: kq = 
{'yf = 2 + i, 7^ = 1 + 2z, 7]^ = 1 — i, 7^ = 1 — 2i). Case c: kq = 1 (tlie anticlockwise direction, 
7^*" = — 1 + i, 74" = — 1 + 2i, 7-f = — 1 — O.li, 7^7 = —1 — 0.2i). Case d: kq = —1 (tlie clockwise 



direction, 7^*" = 1 + i, 7^ = 1 



7i 



1 - 0.3i, 74- = 1 



Due to the fact that A2(— ?/) = Ai(7/), the function e(7/) is odd, and the index of the function l{r]) is 
an even number, ind/(7/) = 2ko. If kq = 0, then the function s(r]) is continuous at the point r/ = 0, 
otherwise it is discontinuous at ry = 0. Notice also that 



Ko 



1 

2^ 



[argAi(r?)]- 



00 
00 



1 

2^ 



arg 



Ai(^) 



A2(^) 



(4.19) 



± 



Numerical computations implemented for different sets of the problem parameters ko, /3, 7^ and 
74 show that Ko may be 0, 1, or -1. Figure 2 shows the graphs of the function l{rj) for < rj < +c« in 
the plane (Re/(r/),Im/(r/)) for some values of the problem parameters in case (i) when both zeros of 
the polynomial 5o(?7) lie in the lower half-plane. The graphs of the function l{r]) for some values of the 
problem parameters in cases (ii) when the two zeros of (5o(??) lie in the opposite half-planes are given 
in Figures 3a and 3b. Figures 3c and 3d illustrate case (iii) (both zeros lie in the upper half-plane) 
when index kq vanishes and equals -1, respectively. In the last case rj traverses the contour around 
the origin in the negative direction. 

Due to the choice of the branch of the logarithmic function log /(ry), if ko = 1, then e(O^) = =F27ri, 

e(0 ) = ib27ri when kq = —1, and e(±0) = in the case kq = 0. Transform now the integral (I4.12P 

into the form 

w /""^ e{t)tdt 



i^o{'n,w) 



27riJo ,/f(t){t^ - r^^y 



(4.20) 
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Figure 3: Parametrically defined function l{i]), < r] < +00, in cases (ii, section 3): a, b and (iii): 
c, d, when ko = 1 + O.li, j3 = 7r/4. Case a: kq = {it ~ ^-^ + 0-5«, 7^ = 1 + i, 7f = 1 — 0.5i, 
7^7 = 1 + i). Case b: kq = 1 (the anticlockwise direction, 7j^ = 1 + z, 7^ = 1 — z, 7-i~ = 1 — i, 
7^ = 1 + i). Case c: kq = (7^^ = 1 + 3i, 74" = 1 + 4i, 7j~ = 1 + z, 7^ = 1 + 2i). Case d: kq = —1 
(the clockwise direction, 7^*" = 1 + 0.5i, 7^'" = 1 + z, 7]~ = 2 + 0.5i, 7^7 = 1 + 2i). 



This integral is bounded as r/ — )■ when kq = and has a logarithmic singularity otherwise. 



il^o{r],w) = {-iy ^KQ\ogr] + 0{l), ?? ^ 0, (r/,u;)Gq, j 
It can be directly verified that the function 

f-qij W + ^ dt 

2^ t~ri 



1,2. 






(4.21) 



(4.22) 



is bounded at the point rj = 0. Here, q^j are the two zero points of the surface, qqj = (0, {—^y^^\/f{0)), 
and qij = {{ — iy~^po, ( — l)-'~^\//(/9o)) £ Cj, where po is an arbitrary fixed point in the upper half- 
plane which coincides with none of the branch points aj {j = 1,2,3,4), the point rjo, and the roots 
of the polynomials di{r]) and 5o{r])6o{—r]). Due to the symmetry of the points qn and qi2 we can 
rewrite the function i{jQ(ri,w) as 



4^oir],w) = 'ipo{'n,w) + Ko 



po 



rjdt 
t"^ — rf' 



+ w 



po 



tdt 



^/m{t'-r]') 



(4.23) 



The function exp{ijjf)(rj,w} satisfies the boundary condition of the RHP (j4.8p . However, at the 
infinite points of the surface it has essential singularities due to the second order pole at infinity of 
the function (j4.23p . In order to remove this singularity, we cut the surface *H along the loops a^ and 
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-a 




Figure 4: The canonical cross-sections a,- and b,-, j = 1, 2, 3. 



hj [j = 1,2,3) which form a system of canonical cross-sections (Fig. 4). The loops a^ intersect the 
loops hj at one point and do not intersect the other loops. The cross-sections ai and a2 and bi and 
b2 are symmetric with respect to the origin. The curves ai, a2 and as are closed contours which 
join the branch points — oi with —02, oi with 02 and 03 with 04, respectively, and lie on both sheets 
of the surface. The loops bi, b2 and h^ join the points —02 with 03, —03 with 02 and —04 with 04, 
respectively, and pass through the infinite points of the surface. The solid lines in Fig. 4 correspond 
to the parts of the loops on the sheet Ci, while the broken lines show the parts of the loops on the 
second sheet. Because of the symmetry of the surface, to remove the essential singularity at infinity, 
it suffices to add to the function ip{ri,w) only the symmetric integrals. 



'po 



PI 



+mo f +mo f +no f +no f J -^. (4.24) 



Here, the point pQ = (ctq, ^/(ctq)) is an arbitrary fixed point in the sheet Ci, while the point 
Pi = i(^i,w{ai)) may lie on either sheet of the surface and has to be determined. The numbers mo 
and n-o are integers to be fixed. The contour poPi is a smooth simple curve on the surface that does 
not cross the loops slj and bj {j = 1,2). 

The function 1^(7], w) admits a simple alteration 



w 



^iv,w) = -— 



e{t)tdt 



27TiJo Vmit^-r]^) 



+ 



Kq 



PO 



rj + 



wt 



dt 



pi 

+mo (t +no 

Po "'32 "'b2 



V ^mjt'-ri' 

w + ^ tdt 



e t^ 



7] 



2' 



(4.25) 



The analysis of the solution at the infinite points of the surface shows that both functions, ip{r], w) 
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and x(^! w) = exp{^p{r], w)}, are bounded as r/ — )• oo if and only if the following condition holds 

1 f°° e(t)tdt fPo tdt fP^ tdt f tdt f tdt , , „ , 



2TTiJo ^/J(T) Jo y^W) Jpo i ^aa C /ba ^ 

4.3 Jacobi inversion problem 

The condition (|4.26|) can be reduced to a genus-1 Jacobi inversion problem solvable by inversion of 
an elliptic integral. First, we express the integrals over the loops a2 and b2 

tdt rl dt 



a2 



e Jal y^Ut)' 



b2 ? Jal ^Mt) Joo ^Mt) 



„2 



tdt r°° dt f'i dt ,^ „, 

+ / ^=== (4.27) 



through elliptic integrals. Here, f{t) = f*{t'^), and the branch ^/*(?7*) is fixed by the condition 
\/JJj]*) '^ 2ir/^sin^/3, r/^, — )• oo, in the r/*-plane (??=k = "rf) cut along the two lines joining the 
branch points, af with a^ and a\ with a\. By making the homographic transformation ()19p t = 
(«! + a2''")(l — I^t)~^ we fix the parameters ai, a2 and /i such that the points af, 02, 03, and of 
are mapped into the points 1, 1/k, — 1/k and —1, respectively {k is to be determined). It directly 
follows that 



t-ai _(l-/i)(l + r) t-af _(l + /i)(l-r) 



al-al 2(1 -/ir) ' a| - af 2(1 - /ir) ' 

and 

t — a| (1 — /i/K)(l + kt) t - 02 (1 + /i/K)(l — kt) 



(4.28) 



2(1 -/ir) ' a2-a2 2(1 - /ir) ' ^^'^^^ 



Equations (|4.28|) imply 

t-al _ (^-1)(1 + t) 



On putting in this equation t = 02, t = 1/k and t = 03, r = —1/k and denoting /^o = (/U — l)/(/i + l), 
/^o = (/^ — l)/('^ + 1) we obtain the system of equations 



/io _ 02 ~ '^4 

2 2' 

kq a2 — af 



Aio^o = -4 1. (4.31) 



It has two sets of solutions 



I {al - al){al - al) ^^ ^ ^ K4-<^i)i4 



/^0-±\/7-2 -2T772 TIT' '^0-±\/772 _2n/_2 3T- l^-^^i 



Each of them defines the parameters /i and k, 

// = - , K=- . (4.33) 

1-^0 1 - KO 

By expressing t from the two relations (j4.28p and adding the formulas obtained we have the homo- 
graphic transformation sought 

^ al + aj jal - al){T - fi) 

2 2(l-^r) ' ^ ' ' 
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We next differentiate the two equations (|4.29p and find 

dt = -^ ,,/ ' ^-^-^ -dr. (4.35) 

2(1 — /iTJ^ 

On multiplying the four relations in ()4.28p and (j4.29p and using ()4.25p and ()4.35p we deduce 

"' .hf „. t ■ ("6) 



where 



1 

isin^/3y {(12 ~ '^3) ('3^1 ~ (^i) ' 



h = —-ni\lTj2 — 12^712 — :2T' (^•2'^) 



e is real, |e| = 1, and v{t) = ^(1 — t'^){1 — k^t'^) is the branch of the two-valued function w^ = (1 — 
T'^)(l — K^T^) in the r-plane cut along the segments [— 1/k, —1] and [1, 1/k] which meets the condition 
u(0) = 1. We wish now to compute the integral over the loop a2 in (|4.27p . Since v{t) = —i\v{T)\ as 
tG [1,1/fi:], from (gSSD, 

tdt , /■!/« dr , , , , 

= eh —— = ehiK', (4.38) 

/aa C Jl v{r) 

where K' = K(\/l — k"^) and K = K(k;) is the complete elliptic integral of the first kind. To fix 
the sign of e, it is sufficient to consider the case when all a^ are real and af < a^ < "^3 < 04. 
The branch \/fJji7) has been fixed by the condition ^/Jjjh) ~ 2ir/^ sin^ /3, ry^, — )• 00, and therefore 
^/fJjh) = «|a/ /*(??*)! as ry* G [oijol]- This implies e = —1. 

To compute the integral over the loop b in (j4.27p . we note that v{t) = —\v{t)\ as 1/k < r < +00 
and derive 

/ ^-^ = -2h r ^ = 2hK. (4.39) 

Denote next 

^i = \^ -, a± = ^—^, j = l,2. 4.40 

fiaj + a_ — a+/i 2 

On making the transformation ()4.34p in the second integral in ()4.26p we derive 

PI tdt _ h f /"O dr /""^i dr \ 

p, T""2 1^4^ + 70 H^J ^ • ' 

in the case pi G Ci and 

fPi tdt h ( /"-i dr /-o dr r^ dr \ , _ , 

(4.42) 



Po 



^ 2 \ y^o ■v(t") J~i v{t) Jo v{t) 



in the case pi G C2. Now the Jacobi inversion problem ()4.26p may be put in a more convenient form. 
Let first pi G Ci. Then 

d + 4noK - 2imoK' , (4.43) 



/o w(t) 
where 



1 /-^ e(t)tdt ^ 2ko r^o trft ^ n J^_ ^^_^4^ 



vri/iJo V7(*) ^ -'o VTW -'o v{t) 
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By inversion of the elliptic integral in (|4.43p we find ai = snd and therefore the affix ai of the point 
pi is given by 

ai = ±Ja, + ^^(^1^. (4.45) 

Because of the symmetry, either sign leads to the same solution. The integers ttiq and no can be 
directly found from (j4.43p 

Im[(J,-j)K] _ Re[(I.-jp] 

° 2Re[KK'] ' ° 4Re[KK'] ' ^ ' ^ 

where /* is an elliptic integral of the first kind 

r"^! dr 



1^=1 ^ = F{sm-\snd),K). (4.47) 

Jo v{t) 

If at least one of the numbers mo, riQ is not integer, then pi £ C2, and no, m,o are both necessarily 
integer. In this case the Jacobi inversion problem (14.26P become 

f^i dr - / 1\ 

-/ — - = (i + 4no + -K- 2imoK'. (4.48) 

Jo v[t) \ 2/ 

The affix of the point pi is the same as in the previous case, while the integers rriQ and no are 

QlIIGrGIlt 

M(/.+j)K| l_Re[(J.+jm 

2Re[KK'] 2 4Re[KK'] ^ ' 

This completes the solution of the Jacobi problem. 

5 Solution to the vector RHP 

5.1 Matrix factorization and its analysis 

As a matter of utility, it will be desirable to have the solution to the RHP (j4.8p on the Riemann 
surface fH expressed in terms of two functions defined on the r/-plane 



where 



X{'n,w) = ey.^{^i{ri) + w^2{ri)}, (5.1) 

2 ?7 + Po 2 r/^ - cr^ 

Because the solution to the Jacobi problem, the point cxi and the integers niQ and no, satisfy the 
condition (j4.26p . it is possible to alter formula (15. 2p for the function V'2(??) as 

1 f 1 /"^ e(t)t^dt fpo t^dt 



v^\2mJo vmit^-v') Jo ^/mit^ - V') 
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which can conveniently be used for large rj. The analysis of this formula shows that the function 
V'2(^) is bounded as r/ — )• oo, and therefore the essential singularity of the solution x(??; w) has been 
removed. 

With the solution to the RHP (I4.25P on the surface D\ at hand, we may write down formulas (I4.9p 
for the factor-matrices X^{rj) and their inverses in terms of functions defined on the r/-plane 



[X±(7?)]-i = e-^tin) cosh(/i/2(,^)v,±(r?))I - -^^^ sinhif/\r^)^f{r^))Q{rj) 



.± 



T] G C=^. (5.4) 



It is essential for the solution of the vector RHP to study the behavior of the matrix X{r]) = X^{7]), 
r] G C^, at the points po = {o'o, \//(o"o)) G Ci and pi = {ai,w{ai)) G *K. For this purpose we use the 
representation ()4.9p . From (jS.ip . ()5.2p . 

X{ri, w) ~ °^ , r/-;>±cro, (r/,u;)G Ci, 

X(??, -■w) ~-E^2> ^-^±cro, (r?, -u;)gC2, (5.5) 



where Er,- are nonzero constants. This yields 



Xir])^^^Y^, V^±ao, (5.6) 

where y^ = y(ib(To, ^/(ibo"o)) is a rank-1 2x2 matrix and E'i are nonzero constants. 
If pi is a point of the first sheet of the surface $H, then 

x{r],w)r^Et^{r]Tai), x{r],-w) ^ Ef^, V ^ ±<^U (r?, (-l)^"^^) G q, j = l,2. (5.7) 

In the case pi G C2, 

xir],w)^Et^, x{ri,-w)^Et^{r,^ai), ?? ^ ±ai, {r,,{-iy-^w) e Cj, j = 1,2. (5.8) 

Here, E^- {s,j = 1, 2) are nonzero constants. Then, regardless of whether the point pi belongs to the 
first sheet or the second one, the inverse matrix [X{ri)]^^ has a pole at the points ibo"i. Assume now 
that ^(7?) is an order-2 vector whose components have certain nonzero limits at the points t] = zbui. 
Then since Y{r],w) is a rank-1 matrix, we have 

[X(r?)]~i^(r/)~-^( i ), 77 ^±^1, Pieq, j = l,2, (5.9) 

where Ef =const, sf = n{±ai)[l{±ai) + wf]~'^, and wf = {-iy-^y/f{±ai), for pi G q, j = 1, 2. 
At the same time, if ^{r]) is an order-2 vector such that 

^(^)~^^(i), ^^±^1, (5.10) 
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where Ej^ are nonzero constants, then 



Ef (I 



where ^(7?) is an order-2 vector bounded as ?? — 7. itiTi. 

We next examine the matrix X{rj) and its inverse as r/ — >■ ipo- If kq = 0, then both functions, 
V'i(??) and V'2(^)) are bounded at the points ipo, and therefore the matrices X{rj) and [X(r/)]~^ are 
bounded as well. If |«;o| = 1) then 

Xiv, w) ^D+, T]^ po, (r?, w) G C2. (5.12) 

At the points with affix —po, 

Xiv, w) r^ D^, T] ^ -po, {7], w) G Ci, 

x{v,w)^D2{7] + po)-'''\ V^-PO, ir],w)£C2. (5.13) 

Here, D- are nonzero constants. At the points rj = itKopO) the matrix X(ri) and its inverse accord- 
ingly behave 

X(7?) = 0(l), [X{r,)]-'r^{rj-KoPo)-'Y^'Hv), V ^ ^^oPo, 

X{r])r^{rj + KoPo)-^Y'^^Hv), [X{rj)]-^ = 0{1), t] ^ -Kopo, (5.14) 

where Y^^'[kqpq) and Y^^\—koPq) are rank-1 matrices. 

Before proceeding with the solution of the vector RHP, we need to factorize the index-0 even 
function y/A(r]). It is given by 

'X{^) = P-M rjeL, (5.15) 

P iv) 

where 

Referring now to the representations (|3.20p and (j4.7p . we split the matrix G{r]) as 

G(r?) = ^i^p+(r?)X+(r?)[p-(,?)X-(r/)]-\ r? G L. (5.17) 

We remind that 60(7]) and G^{ri) are the polynomial and polynomial matrix given by (j3.14p and 
(|3.17p . respectively. 

5.2 Vectors $+(77) and ^"(r/) 

To determine the vectors $+(77) and <I>^(r/), we substitute the splitting ()5.17p into the boundary 
condition (|3.8p of the vector RHP and use the formulas 



On factorizing the polynomial ^1(77), we have 

^o(r?)[X+(r?)]-iGi(77)cD+(r?) _ (r? - ti)(r? - ta) 
5^a2{ri + ti){r] + t2))p^{r]) P'iv) 
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[X-(r/)]-icl>-(7?), 77 GL. (5.19) 



Here, itti and ±t2 are the four zeros of the polynomial di{r]) = Q2(??^ — ii)(^^ ~ ^Dj ^^^ Imtj > 0, 
j = l,2. 

We begin with case (i) when the two zeros of the quadratic polynomial So{r]) lie in the lower 
half-plane, So{r]) = —{r] + Ti){r] + T2)sm'^/3, Imr^ > 0, j = 1,2. On applying the principle of analytic 
continuation and the Liouville theorem, we have 

.(. + r.)(..r.)[xn.,)]-'6'W^n.) _ O-'.K--'^) ,^-,,,)]-.-,,,) = ..(,,), , e c. 

(5.20) 
where v = —5^ csc^ /3, and i?(r/) is a rational vector-function. Note that due to (j5.9p the vector i2(r7) 
has simple poles at the points itcri. Also, because of ()5.14p the vector R{i]) has simple poles at the 
points rj = kqpq in the case kq = ±1 and is bounded if kq = 0. In addition, it has to have simple poles 
(the geometric optics poles) at the points rj = ±770- Since the vectors $ (ry) vanish, the matrices 
[X^(r/)]~^ are bounded at infinity, and the elements of the matrix G^(ry) are degree-2 polynomials, 
the vector R{r]) has a pole at the infinite point of multiplicity 1 if kq = and multiplicity 2 if 
kq = =tl- The most general form of the vector R^r]) in the former case is given by 



PJV) 
iv^ - (^i){v^ - Vo)'' 



Piv) = ,^, _2JL 2T' (5-21) 



where P{'i]) is an order-2 vector whose components, Piit]) and P2{v)j ^^^ degree-5 polynomials. In 
the case kq = ±1, 

P{ri) 

where the components Pi{t]) and P2{v) of the vector P{ri) are degree-6 polynomials. The vectors 
$^(r/) can be deduced from (j5.20p 

^+(„) = ^pHri)GHn)x+{ii)R{v) + 

$-(.) = ^f^^^^^T^, .€C-. (5.23) 

The solution has twelve arbitrary constants, the coefficients of the degree-5 polynomials Pi{r}) and 
P2{'n) ™- the case kq = and fourteen constants if kq = ±1. 

In case (ii), when one of the zeros of the polynomial ^oirj)^ ri, lies in the upper half-plane and 
the other, T2 G C2, (^o(^) = —{v — Ti){'n + ^2) sin^ /3, the boundary condition of the RHP gives 

K-, + r,)[xn,r'6'W*-W _ (»-'.)C-*=)[x-(„r..f-(,) = i,M. „<,c. (5.24) 



iv + ti)ir] + t2)p+{v) {v-n)p iv) 

The rational vector R{r]) is given by (j5.2ip . where Pi{rj) and P2{'n) ^'^^ polynomials of the fourth 
degree in the case % = and by (j5.22p with the fifth degree polynomials in the case kq = ±1. The 
rearrangement of the factors of the polynomial (^o(^) produces the new solution 

{ri + T2){r] - ti){r] - t2) ' 

^-( h^n)p_(^ix:i5m, ,,c-, (5.25) 

(j] - ti){r] - t2) 
which has either ten, or twelve arbitrary constants depending whether kq = 0, or kq = il- 
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In the third case, when both zeros of the polynomial (5o(fy) lie in the upper half-plane, Sq^tj) = 
— {rj — Ti){ri — T2) sin^ /3, we have 

{tl + h){ri + h)p+{ri) (>)-n)(Jf-T2)p-(t,)' *"' '"' ' * ' 

As before, the rational vector R{r}) has the form (j5.2ip for kq = and (15.221) for kq = ±1. However, 
the components of the vector P{r]) are degree-3 polynomials if kq = and degree-4 polynomials if 
Ko = ±1. The solution has eight arbitrary constants in the former case and ten constants in the 
second case. It may be written as 

^~^r,) = ^^--^)^^---)p-^^)^-^^)R^^) ^ ,eC-. (5.27) 

^' {v-tiKr]-t2) ^ ' 

5.3 Symmetry conditions 

In general, the solution derived in the previous section does not meet the symmetry condition 
$"'"(r/) = $~(— ry), rj £ C We begin with case (i). In order to satisfy the symmetry condition, 
we rewrite it as 

fi(-„)= :'";">'" ^''^''f;,""', , ix-(-„r^a\,)xu,m:„). v<=c-. (5^28) 

iv + n)(r? + T2)[ri - ti)[r] - ta) 

Here, we used the relation p^{t]) = [p~(— ??)]~^, ^ S C^. We wish now to simplify the matrix 
U{r]) = [X~ {—r])]~^G^ {r])X'^ (rj). It will be convenient to represent the function V'i(^) as the sum of 
the even and odd functions 



— log— , V'lei?? = - log ^ 2 



ipiiv) =i^ioir])+ipie{v), ^io(^) = — log;:— — , V'le (??) = - log 35 ^, (5.29) 



(the function V'io(??) = if kq = 0) and use the following notations: 



c = cosh[^J f{±r])^lj2i±v)], s = sinh[y^/(±??)V'2(±??)], 
By referring now to (15. 4p . we transform the matrix U{r]) as 

^( ) ^^2,.(.) f ^-^^Yyj ---///') (all s\A f ^+^y^ ''"''/ ffj] ■ (5-31) 

\-sn^/Vf c + sl-Hf ) \92i 922 )\ sn+I^Jf c-sl+/VfJ 
On using next the directly verified identities 

{1+ - l-)9n + n+9i2 - ^-oli = Oi "i+S'ii - {h + l-)9i2 - ^-922 = 0' 

-n-9\i + [1+ + 1^)921 + "'+5'22 = 0> -n~9i2 + "1+521 " (^+ " ^-)5'22 = 0> 

c^ah + -j{-l-l+9ii - l-n+9l2 - l+m-gli - m_n+5r22) = 9ii, 
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^2 

c^5i2 + y (-^-"T'+S'ii + 1-1+912 - m-m+gli + l+m-gl2) = 5i2; 

s2 
c^92i + -j{-hn-g\i - n_n+g\2 + l^l+g\i + l^n+gl^) = gli, 

c'^922 + -ji-'m+n-gli + l+n^gl2 + l-m+gl;^ - Z-/+5'22) = 522' (5-32) 

we ultimately deduce that U{r]) is a polynomial matrix and it coincides with G^{r]) that is everywhere 
in the upper half- plane 

[X-{-r^)]-'G\v)X^{r^) = fH^Y cHr^). (5.33) 

\V + PoJ 

We show next that [/o'^(^)]^ is a rational function. Consider the integral 

where A(t) is defined in (j3.25p . On factorizing the function A(r/) we obtain 

Hence we have to deal with the integral 

which can be evaluated explicitly 

Thus, we have proved that [p^{t])]'^ is a rational function, and it has the form 

(?? + il)(?? + t2) 

For the present purpose, we transform the symmetry condition (j5.28p as 

R{-rj) = J^ fH^fRY" G\t])R{t]) (5.39) 

or, equivalently, 

= (-')--g;C)^("). (5.40) 

where the components of the vector P{r]), Pi{rf) and ^2(^)5 are polynomials 

|reol+5 |ko|+5 

j=0 jr=0 

whose coefficients are to be determined. On replacing r] by —77 we reduce the condition (|5.40p to 

22 



which is equivalent to the condition (|5.4U|) since the matrix G^{r]) possesses the following property: 

G\ri)G\-rj) = di{r])I. (5.43) 

The symmetry condition (j5.40p asserts that the coefficients of the polynomials Pi (77) and P2 {rj) can 
be determined from the two equations 

Qi{r])=0, Q2{r])=0, (5.44) 

where 

Qjiv) = sin2/3(r? - ti)(r? - t2)P,{-v) + {-irS4g],{r,)P,{i^) + g]2{v)P2{v)], 3 = 1,2, (5.45) 

are polynomials of degree 6 + |«;o| (the terms Ty^+I^'ol in both polynomials have zero-coefficients). 
Consider first the case kq = 0. To define the coefficients, we compute the derivatives Qj {ff}, 

m = 0, l,...,6;j = 1,2. The first six equations Q\ (??) = {m = 0,1,..., 5) express the 
coefficients h^-m through the coefficients of the polynomial Pi (r/) 

b2j = vo{-a2jvi + a2j+ii^2)^ 

b2j+i = i^oi-2a2jsm^ /3 -a2j+ii^i +a2j+2i^2)^ J = 0,1,2, (5.46) 



where 



'' = 2cosPUy-V'^ + k^sin^^y -f = 7r-7r±(ti+t.)sm^/3. 



zy± = 7i 74 + k'^ cos^ /3 ± tit2 sin^ ^, uq = 0. (5.47) 

The equation Qi{rj) = yields 

ao [27f Tzf + k"^ cos^ l3 - tit2 sin^ /3] = 0. (5.48) 

It is directly verified that 27^7^ + k"^ cos^ /3 — tit2 sin^ /3 = and therefore oq is free. The other 
seven equations Q2 = (j = 0, 1, . . . , 6) are also identically satisfied provided the coefficients bj 
{j = 0,1, . . . , 5) are chosen as in (15.46p . and ti, ^2 are the two zeros in the upper-half plane of the 
even degree-4 polynomial di{7]). 

A similar result holds for the case |«;o| = 1- The equations Ql ^ {r]) = {m = 0,1,..., 6) 
express the coefficients be-m through the coefficients of the polynomial Pi{f]) 

&2j = z^o(-2a2j-isin^/3-a2jZ^+ + a2j+iJ^2')' i = 0,1,2,3, a_i = 07 = 0, 

62J+1 = fo(-a2j+iz^r + a2j+2i^2^), j = 0, 1, 2. (5.49) 

In the case kq = ±1, we have 2^^^^ + k"^ cos^ /3 + tit2 sin^ /3 = 0. Therefore, the equation Qi{rj) = 0, 
which is equivalent to 

ao[27r74: + k"^ cos^ /3 + ^1*2 sin^ /3] = 0, (5.50) 

is satisfied for any gq. If the coefficients bj are defined by (I5.49p . then Q2{'n) = 0- 

Consequently, in case (i) , the solution to the vector RHP possesses 6 + | kq | arbitrary constants 

aj (j = 0, 1, . . . , 5 + |ko|) and satisfies the symmetry condition ^'^{■q) = $~(— r/). 

Consider next case (ii) when the two zeros of the polynomial (5o(^) he in the opposite half-planes. 

Since the relation (I5.33P is invariant with respect to the location of the zeros of the polynomial (5o(f?), 

we employ it again and deduce from (j5.25p 

Ri-V) = -7 tZ TT ("^y G\v)R{v)- (5.51) 

(77 -ti) (77-^2) \V + PoJ 
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The rational vector R{r]) is given by (|5.21|) if kq = and by (|5.22|) if kq = =tl) and 

|ko|+4 |ko|+4 

Piiv) = E «J-^'' ^2(r/) = E ^jV- (5.52) 

j=0 j=0 

The symmetry condition (3.89.3) reads for the polynomials Pi{ri) and P2{'n) 

On following the procedure described for case (i) we obtain for kq = 

b2j = i'Q{-2a2j-i SVC? P - a2jut + (^2j+iJ^2)^ J = 0' 1' 2, a_i = 05 = 0, 

hj+i = J^Q{-a2j+ii^i + a2j+2i^2)^ J = 0, 1- (5.54) 

If |ko| = 1, then the coefficients 62J and ^2^+1 coincide with the corresponding coefficients in case (i), 
K = 0, and are defined by (jSliGj) . In case (ii), the solution is given by (jOSj) . ([OT]) . ([02]) . (f532|) . 
(|5.54p and ()5.46p . It satisfies the symmetry conditions and possesses 5 + |ko| arbitrary constants 
ao,ai, . . . ,a4+|Kg|. 

Finally, in case (iii), when the two zeros of the polynomial 5o(^) lis in the upper half-plane, the 
symmetry relation coincides with the one derived in case (i) and is given by (|5.39p . The polynomials 
Pi (7/) and P2{'n) in the representation of the solution (I5.20p have the form 

|ko|+3 |k()|+3 

^i(^) = E «^-^'' ^2(r?) = E ^^■^'' (5-55) 

i=o j=o 

where the constants aj (j = 0, 1, . . . , 3 + |ko|) are free and if kq = 0, 

&2j = VQ{-a2jh'i + a2j+ii^^), 

62J+1 = VQ{-2a2j sin^ /3 - a2j+iVi + a2j+2J^2), j = 0, 1, 04 = 0. (5.56) 

If I Ko I = 1 , the coefficients are defined by ()5.54p . The solution derived has 4 + 1 kq I arbitrary constants 
ao,ai, . . . ,a3+|K(,|. 

5.4 Additional conditions 

Since the vector polynomial P{r]) meets the condition (j5.40p in the cases (i) and (iii) and the condition 
(I5.53|) in the case (ii) for all rj, and the left-hand side in (j5.40p and (15.53P is bounded at the points 
ti and t2, due to ()5.23p . ()5.25p and ()5.27p the vector $+(77), regardless of which case is considered, 
(i), (ii), or (iii), has removable singularities at the points ti and ^2 (Imtj > 0). 

Now, the property ()5.9p of the matrix [X(77)]~^ asserts that at the point ai the polynomials Pi{r]) 
and P2ir]) satisfy the relation 

stPi{c7i) = P2{<Ji). (5.57) 

Because of the symmetry condition this relation is necessary and sufficient for the removal of the 
simple pole of the solution at the points zto"!. 

According to the property ()5.6p the solution to the vector RHP has simple poles at the points 
ibiTo. Since (Tq is the affix of an arbitrary fixed point po on the first sheet of the surface *H we 
may assume that IniCTo > 0. Since Y'^ in (|5.6p is a rank-1 matrix, by the symmetry property of 
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the solution, the points zbuo are removable points of the vectors ^^{rj) if and only if the following 
condition holds 

Y+Pi{ao) + Y+P2{ao) = 0, (5.58) 



where {¥{^,¥{2) is the first row of the matrix ¥{ao, y/fjoo)). 

The solution derived has geometric optics poles at the points itr/o (^0 = k sin 9o sin /3) with the 
residues of the vector <&+(r/) given by ()3.7p . This yields two extra equations for the determination 
of the six coefficients Oj . They are 



where 



ESi = iEo, ES2 = iZHo, (5.59) 



2%(^o - f^i)(% + ri)"i(r?o + T2)"2(% - ii)(% - ^2) ' 

Sj = [9]i{m)Xu{m) + 9hivo)X2i{m)]Pi{vo) + [g]i{m)Xi2{vo) + g]2ivo)X22{m)]P2{r]o), (5.60) 

where X{r]) = (Xm„(r/)), m,n = 1,2, and ai = 02 = 1 in case (i), ai = 0, 02 = 1 in case (ii), and 
ai = 02 = in case (iii). 

Let Ko = 0. There are s arbitrary constants oq, ai, . . . , Og-i, where s is equal to 6, 5 and 4 in the 
cases (i), (ii) and (iii), respectively. At the same time, there are four additional conditions regardless 
of which case is considered. Thus, the solution has two arbitrary constants in case (i), one constant 
in case (ii) and it is unique in case(iii). 

When Ko = 1, according to (I5.14P and (I5.22p . the vector ^^(rj) has inadmissible pole at the point 
r] = pq £ C"*" (the point po was chosen to be in the upper half-plane) due to the pole of the vector 
R(r]). The vector $~(r/) on the other hand has a simple pole at the point rj = —po. Since y_ in (j5.14p 
is a rank-1 matrix, for the point r] = —po S C~ being a removable point of ^^(^) it is necessary and 
sufficient that the following condition holds: 

YuPii-po) + Yi2P2{-po) = 0, (5.61) 

where Yu and Y12 are the elements of the first row of the matrix Y^^'{—pq). If the condition (I5.6ip is 
satisfied, because of the symmetry condition ^~^(r]) = $"(—77), the point rj = po G C"*" is a removable 
point of the vector $+(77). 

In the case kq = —1, the matrix X{r]) has a simple pole at the point rj = po, and the vector R(t]) 
has a simple pole at r; = —po- Therefore, the solution $ (ry) has simple poles at the points ±po- 
Similarly, to the previous case, it can be removed by the condition 

^ii^i(po) + ^i2^2(po) = 0, (5.62) 

where Yu and Y12 are the elements of the first row of the matrix Y'^^'Ipq). 

Because in the case \kq\ = 1 we have s + 1 free constants, ao,ai, . . . , fflsi the number of free 
constants exceeds the number of additional conditions (it is five) by two in case (i) and by one in 
case (ii). In case (iii), the number of free constants and additional conditions is the same. 

This completes the solution procedure for the RHP 1. Since the structure of the matrix coefficient 
G{r]) for the RHP 2 is the same as in the RHP 1 (the five parameters of the problem, 7^ , 74 and 
P, need to be changed according to the transformation ()2.19p ). the same method is applied to the 
RHP 2. As in the case of the RHP 1, the solvability of the RHP 2 is governed by the polynomial 
So{r]) which now becomes 

~6oiv) = {rf - kl) cos2 /3 - (7? - 7+)(r? - 7+). (5.63) 
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To summarize the result proved, we introduce a class of solutions. We say that 4'j{x, y) G ^(VFg) 
if <^,{x, y) e C\We) n C\We) (i = 1, 2) and 

\<l>^-Ei-El-El\<K^e-'P, 
\Z-^(j)2-Hl-Hl-Hl\<K2e~^P, < p < oo, 0<6'<|, (5.64) 

where e is a small positive number, Ki and K2 are finite positive constants, {El, HI) is the incident 
wave ()2.2p . {El, HI) is the reflected wave, and {El, HI) is the surface wave. 

Theorem 5.1. Let a right-angled concave wedge he illuminated by the incident wave i2.1\) and the 
material properties of the wedge faces be described by the impedance boundary conditions (12. 6|) . Let 
a = 'K 12 and 72 =73 = 0. Then, in the class of functions A{We), the boundary value problem (12.8)) . 
( 12. 9|) has n + 1 solutions, where n is the number of zeros in the lower half-plane of the characteristic 
degree-4 polynomial 

'5o(^)^o(r?) = [{V^ - kl) cos2 /3 - (7? - 7r)(r/ - ^^)][{ri' - kl) cos^ /3 - (r? - 7i+)(r/ - 74+)], (5-65) 

where feg = fesin/3 and Im/c > 0. In particular, the problem has a unique solution if and only if all 
the four zeros lie in the upper half-plane. 

5.5 Scalar case: (3 = it/2 

As a test for the solvability results of Section 5.4, consider the scalar case when a = /3 = 7r/2. The 
matrix coefficient G{i]) of the RHP 1 becomes the diagonal matrix 

G{i]) = diag{(7f + r]){'y^ - r/)-\ (7^ + rj){-f^ - r/)"^}, (5.66) 



and the vector RHP (j3.8p reduces to the following two scalar RHPs: 

Find two pairs of functions ^t{r]) and $7(r/) analytic in C"*" and C~ , respectively, except for the 
geometric optics poles r/ = itr/o with residues defined by (13.7]) and continuous up the real axis L. On 
the contour L, their limit values satisfy the boundary conditions 

<^+{rj) = l2^±l<l>j{n), r,eL, j = 1,2, (5.67) 

■' 7j_ — T] ■' 

where 7i_ = 7J~, 72- = 747. At infinity, ^j{i]) = 0{r]^^), and everywhere in the plane ^t{rj) = 

$7 (-??). 

The coefficients of the problems are rational functions, and the solution can easily be derived. 
The polynomial 5o{ri) introduced in (|3.14|) has the form 5o(^) = ~{v ~ li)iv ~ Ia)- ^s before, we 
consider cases (i), (ii) and (iii). If both zeros of the polynomial (5o('?) he in the lower half-plane, that 
is if Im7j_ < 0, J = 1,2 (case i), then the solution to each RHP (|5.67p has one arbitrary constant 
and has the form 

where 

Coi = -Cii?7o + 2i-E'or?o(7r " %)' ^"02 = -CuVq + 2iZHo'no{'y^ - %), (5.69) 

and the constants Cn and C12 are free. 

In case (iii), when both zeros of the polynomial lie in the upper half-plane, Im7j_ > 0, j = 1, 2, 
the solution is unique and has the form 

*+(„) = e%±2), ^jW- Qfa--"' , , = 1.2, (5.70) 
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where 

Ci = ^^^^°^, C2 = ^^^^^^^^. (5.71) 

Vo + 7i Vo + 74 

Finally, in case (ii), when one of the zeros say, 'j^ lies in the upper half-plane and the second 
one 7^7 G C~ , the RHP for the electric field has the unique solution ()5.70p . (|5.7ip (j = 1), while 
the solution to the RHP for the magnetic field is given by (|5.68p . (|5.69p (j = 2) and possesses one 
arbitrary constant. Thus, in each case the results for the test scalar problem are consistent with 
those obtained for the vector case in Section 5.4. 

The RHP 2 reduces to the two scalar problems ()5.67p . where 7j_ needs to be replaced by 7j+, 
and the solution is unique if and only if both Im ^^ and Im 7^ are positive. Otherwise, if Im 7^^ < 
and Im74' < 0, then the solution has two arbitrary constants, and if Im7J^ Im7J^ < 0, then there is 
one free constant. 

Notice that this result is consistent with the uniqueness theorem proved in (llSp for the scalar 
case (/3 = '/r/2) for a wedge W : \6\ < a, a £ (0,7r) by applying the Green theorem to the integral 

(A + k^u)udW (5.72) 

w 

over a curvilinear rectangle W : po < p < pi, —a < 9 < a {0 < po < pi < 00). Here, u is the solution 
to the scalar impedance boundary value problem for the Helmholtz equation in the wedge W. 

6 Sommerfeld integral representation 

In this section, we consider the case when all the four zeros of the characteristic polynomial So{s)6o{s) 
lie in the upper half-plane, and both RHPs have a unique solution. The electric and magnetic field 
can be given in terms of the Sommerfeld integrals as 

Ez{p,0) \ _ Jl_ f ^-ikopcoss I Se{s + 9) 



ZH,{p,9)J -2^i7r' '°'""l,sl(s + e)j'^"' ^^-^^ 

where T is the Sommerfeld double loop with the asymptotes s = tt/2 and s = —3tt/2 for the upper 
loop 7+ lying in the upper half-plane. The starting point of the contour is chosen to be s = 7r/2 + ioo. 
The second loop T- is symmetric to 7+ with respect to the origin. 

6.1 Spectra S'e(s) and Sh{s) 

The purpose of this section is to derive expressions for the functions Se{s) and Sh{s). The functions 
Se{s) and Sh{s) are analytic everywhere in the strip | Res| < 7r/2 + e (e > and small) except for 
the simple poles s = ±0o, and 

res Seis) = Eo, res Shis) = ZHq. (6.2) 

s=6»o s=eo 

The knowledge of these functions enables us to determine the reflected, surface and diffracted flelds. 
We begin with the function Se{s). On applying the inverse Maliuzhinets transform ()20p . 

Se{9 + s)-Se{9-s) = ikosms e'''°P''°'''E,{p,9)dp, (6.3) 

Jo 

and fixing first ^ = and then 9 = tt/2, we obtain the following two equations: 

Se{s) - Sei-s) = ikosins / e^'=°^™""£^(x, 0)(ix, 

JO 
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Se(^^ + s^-Se(^^-s^=ikosmsJ^"e"'^y'''''E,{0,y)dy, (6.4) 

where the functions Ez{x,0) and Ez{0,y) are the sums of the associated functions generated by 
the waves El'^ (case 1) and -E*~ (case 2). In the first case, the integrals in (|6.4p with respect to x 
and y are the functions (/>i(A;o cos s, 0) and 0i(O, ko cos s), respectively, and can be expressed through 
the solution to the RHP 1, the functions ^^ [t]) and $2 iv)- ^^ the second case, these integrals 
are the functions (/)i(/co coss,0) and (/>i(0, /cq coss), respectively, and representable in terms of the 
solution to the RHP 2, the functions "^i {r]) and ^2 iv)- O^ employing the first equation in (j3.4p we 
determine cj)i{kQ cos s , 0) in case 1 and and 0i(O, fco coss) in case 2 (the parameters in ()3.4p have to 
be transformed according to (j2.19p ). In both cases we take r] as ko sin s (because of the symmetry of 



the equation both solutions rj = itfco sin s of the equation ikQ cos s = — y ??^ — kQ produce the same 
result). Equations ()6.4p now have the form 

Se{s)-Se{-s)=Fie{s), Se (^^ + S^ - S^ (^^ - S^ = F2e{s), (6.5) 

where 

Fie{s) = ikfj sin s'^^ [kQ COS s) {jf + /cq coss)[<I>J|'"(A;o sins) — <^|f (fco sins)] 

+-fco COS /3 sins [$2' (^0 sins) + <I>2^(fcosins)], 

F2e{s) = ikQ sins<I'^(A;o coss) — -(7j~ + A;ocoss)[^J^(A;o sins) — ^j~(A;o sins)] 

— -fco cos /3 sin s[^2^(/co sins) + ^^(A;o sins)]. (6.6) 

We will adopt the Fourier transform in the form (I20p and use the following notations for the 
transforms and their inverses: 

1 /'ZOO pioo 

■^TT J —ioo J —ioo 

Fje{t) = -— / F,e(s)e^^*ds, F,e(s) = / F,e{t)e-''' dt. (6.7) 

^^ J —ioo J —ioo 

When applied to the functions 5e(vr/2 it s) this transformation yields 

1 rioo /„- \ pTi-iTt/2 /-n /2+ioo 

-— sA-±s] e'^'ds = —— / SAsV^'^'ds 



•^Vr J —ioo \ ^ / -^TT J -K jl-ioo 

= 5e(±t)e^*"*/2 _ iSQe±i(eo— /2)i, (6.8) 

Here, we used the fact that the function 5e(s) is analytic everywhere in the strip | Res| < 7r/2 + e 
apart from a simple pole at the point s = 9q with the residue given by (j6.2p . From these we find 
that corresponding to ()6.5p we have the pair of equations 

Se{t) - Se{-t) = Fie{t), 

^^(i)e--*/2 _ ^^(_t)e-*/2 = ^2^(i) + 2Eo sin {^ - 9^ t. (6.9) 

Because of the symmetry of these equations, it suffices to determine one of the functions Se{t), Se{—t) 

^ ^ 2isin7rt/2 sin7rt/2 ^ ^ 
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It only remains to apply the inverse transform (|6.7|) and use the integral 



2tano-, <Rea<-, (6.11) 

sin7rt/2 '2 2' ^ ^ 

to obtain an integral representation for the function Se{s). We find 

Seis) = 5(i)(s) + SP{s) + Eo[cot{s - 0o) - cot(s + 9o)], (6.12) 

where 

5(^)(s) = -— ^ / Fie(cr)cot(cr - s)(icr, < Res < vr, 

Si^Hs) = --^ r°° F2e(a)tan(a - s)(ia, -J < Res < J. (6.13) 

27rz J_ioo 2 2 

The function 5e(s) is analytic everywhere in the strip < Res < tt/2 except for the point s = ^O; 
where it has a simple pole with the residue equal to Eq. In order to recover the reflected, surface 
and diffracted waves, we analytically continue the function Se{s) from the strip < Res < 7r/2 to 
the right and to the left to the strips 7r/2 < Res < 27r and — 37r/2 < Res < 0. We have found that 
the solutions to the RHPs 1 and 2, the functions ^j{ri) and ^^(j?), have simple poles at the points 
ibr/o and ±f/o, respectively. It will be convenient to represent them as 



^fiv) = ^Z^, ^f{v) = ^^, J = 1,2, (6.14) 



where $ • (rj) and ^ • (r/) are bounded and nonzero at the points zbryo and ±f/o, respectively. Starting 

with the function Se (s) given by (j6.13p and analytic in the strip — 7r/2 < Re s < 7r/2 we continue 
it to the strip tt/2 < Res < 37r/2 by shifting the contour Ims = to the right by vr. In the strip 
< Res < vr, the function F2e{s) has two geometric optics poles, s = 7r/2 — ^o and s = 7r/2 + ^o- 
In addition, it may have some poles which generate surface waves. Indeed, when < Re s < vr 
and Im(feocoss) < 0, the function ^f{kocoss) needs to be continued analytically from the domain 
< Res < TT and Im(A;ocoss) > 0, where it is analytic (except for the geometric optics poles). The 
same is applied to the functions ^^(fcosins) and ^"(/cQsins) {j = 1,2), when Im(/cosins) < and 
Im(A;osins) > 0, respectively. From the RHPs 1 and 2 boundary conditions we derive 

^ W = :T^-^ ^ W, ^ W = 7-- * [V), f? G C , (6.15) 

where G^irj) and G^{i]) are obtained from G^{r]) and G^{r]), respectively, by the transformation 
(|2.19|) . Here, C{ko sin s) = —iko cos s and ({ko cos s) = —ikQ sin s, and G^ (r/) and & (r/) are polynomial 
matrices. Notice that 

d(A:ocoss) = — (i^(A;osins)(i^(— /cq coss), d{kQsms) = —d^{kQCOss)d^{—kQsms), 

d{k(jCoss) = — d~{ko sin s)d'^{—kQ cos s), d{k()sms) = — d~ [ko cos s)d~^{—ko sins), (6.16) 

where 

d^{t) = f sin^ ^ + (7± + 7±)t + 7f 7± + kl cos^ /3. (6.17) 

Notice that d~{—t) and d'^{—t) coincide with the polynomials (5o(t) and 5o{t) given by ()3.14p and 
(I5.63p . respectively. Introduce next the numbers 

,± _ -{i t+ii) + v^ ^± ^ -{it+ii)-v ^ 

2fcosin2/3 ' ^ 2/cosin2/3 



*1 - n;_ „:. 2 o ' *2 " ^1 T" 5-3 J (6-18) 
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where 

S± = (it + 74^)' - 4 sin^ /3(7f 7± + k^ cos^ /3). (6.19) 

In the case of the RHP 1, we need the zeros of the functions d{kQsms) which are defined by the 
equations cos s = tf2 and sin s = — 15~2 and given by 



s = ±i log(t+2 ± ^(^^2)^ - 1) + 2vrfc', s = 7r/2 T k/2 + i logi-it^^ + V ^ " (^m)^)] + 27rA:", (6.20) 
Also, we require the zeros of the function d{kQ cos s) which are 



s = ±ilog{-t{.2±\/{ti,2)^ - l) + 27rm', s = 7r/2=F[vr/2 + ilog(it]^2 + V ^ ~ (^tz)^)] +2^"i"- (6-21) 

Here, k', k" m' and m" are integers. When the RHP 2 solutions are analytically continued, we need 
to know the zeros of the functions d{kosms) and d{ko cos s). The zeros of the former function are 
given by 



s = ±i log(ti;2 ± ^it^,2? - 1) + 2^A:', s = 7r/2T [vr/2 + i log(-it+2 + V 1 " (^^2)')] + 2vrA;", (6.22) 
while the zeros of the function d{kQ cos s) = —d'^{— cos s)d~ (sin s) are 



s = ±i log(-t^2 V (^m)^ - 1) + 2vrm', s = 7r/2 =F [7r/2 + ^ log(iti;2 + V 1 - (.tl,2)^)] + 2vrm". (6.23) 

Denote Sg • (j = 1, . . . , 722) the poles of the function d{ko cos s) lying in the domain < Re s < vr, 
Im(A;ocoss) < and the poles of the function d{kQsms) in the strip < Res < vr (call these poles 
the surface wave poles). Let 0j (j = 1, . . . , 4) be the geometric optics poles in the strip < Re s < vr, 
the points 7r/2 — 9q and tt/2 + 9q. We can write 

1 rn+ioo ^ 

Sf'{s) = —- — : / F2e{(j) tan(o" — s)da + ^J res [F2e{<j) tan(cr — s)] 

+ E _res^)[^2e(^)tan(a-s)], ^<Res<^. (6.24) 

It will be convenient to use the following notations 

Mt = '^^f°.''"^°^ + ,,2 ■ ,, {(7r ± kosm9o)[i!tikocos9o) - ^^{kocos9o)] 
2kq sm 9o zKq sm 29q ^ 

+ko cos/3cos0o[^2^(^ocos0o) + ^2^(^0 cos^o)] f 7 

N^f^ = res F2eis). (6.25) 



On employing now formulas (j6.6p and ()6.14p we arrive at the analytic continuation into the strip 

(2) 
tt/2 < Res < 3it/2 of the function Se (s). It has the form 

SP{s) = -— ^ / F2e(cr + 7r)tan(cr - s)da + M^" cot(s - 6^0) 

ZTTZ J —ioo 

- M+ cot(s + ^o) - E ^if^ tan(s - sg^^^ ) , | < Re s < ^^ . (6.26) 
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(2) 

To describe all surface waves, we need to continue the function Se further to the right up to the 
line Inis = 27r. In the strip vr < Res < 3tt/2, the function -F2e(s) has one geometric optics pole 
37r/2 — ^0 and poles Sg- (j = 722 + 1, . . . ,722) which may generate surface waves. By shifting the 
contour to the right by tt/2, similarly to (j6.26p we obtain 

S(2)(s) = —\ F2J a + — j cot(CT - s)da + M" cot(s - ^o) 

+ (M- - M+) cot(s + 0o) - E ^?^^ tan(s - ^f ^), — < Re s < 27r. (6.27) 

Analyze next the function S\ (s) given by (j6.13p . It is analytic in the strip < Res < vr and can 
be continued analytically into the strip tt < Re s < 27r by shifting the contour Re s = to the right 
by vr. In the strip < Res < vr, the function Fie(s) has simple poles at the points ^o and vr — ^o- 
In addition, it may have surface wave poles at the zeros of dik^ cos s) in the domain < Re s < vr, 
Im(A;ocoss) < and the zeros of d(fcosins) in the strip < Res < vr. Denote these poles by Sg- 

{j = l,...,ni). The analytic continuation of the function S\ (s) into the strip vr < Res < 2vr 
eventually gives 

Sr{s) = -75-^ / Fie(cr + vr) cot((T - s)da + A+ cot(s - 6*0) 
- A- cot(s + ^o) - ^ ivi]^^ cot(s - s^J^^), vr < Re s < 2vr, (6.28) 



i=i 



where 



i'^i {ko cos ^0 



± = ^-^1 1^0 cos .0; ^ ^ |( + ± kocosOomtikosmeo) - ^HkosinOo)] 

2kQ cos 0Q zKq sm 20q ^ 

— A;ocos/3sin0o[^2^(^o sin^o) + i2^(^o sin^o)] r 1 

iv(J«)= res Fie(a). (6.29) 

After combining formulas (|6.12p , (j6.13p , (|6.26p to (|6.28p we obtain the representations of the function 
Se{s) in the strip < Res < 2vr. Let 

1 rioo 

Ije{s, a) = : / Fje{a + a) cot{a — s)da, a < Re s < vr + a, 

1 rioo 

/jg(s,a) = -— ^ / Fje{a + a)tan{a - s)da, -7r/2 + a < Res < vr/2 + a. (6.30) 

ZTTl J — ioo 

Then we may write 

Se{s) = 11^{s,0) + I^2e{s,0) + Eo[cot{s - 9o) - cot{s + 60)], 0<Res< ^, 

Se{s) = Ilis, 0) + 4(s, vr) + (^0 + M-) cot(s - ^o) 
n'2' 
-(£^o + M+)cot(s + 6'o)-X]^'ir^tan(s-sg^^), |<Res<vr, (6.31) 

j=i ^ 
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where Sg • (j = 1, 2, . . . , 722) are the surface wave poles of F2e{s) in the domain < Re s < 7r/2, 
Seis) = Il^[s, vr) + /*g(s, vr) + (^0 + M" + A+) cot(s - ^0) - (-^o + M+ + A") cot(s + ^o) 

- i: ATij^) cot(. - .(;.«)) - E ivif ) tan(. - sf\ . < Re. < |, (6.32) 

where Sg- (j = 1,2, .. . ,n']^) are the surface wave poles of -Fie(s) in the domain < Res < 7r/2. 
Finally, as 37r/2 < Res < 27r, 

5e(s) = 7fg(s, vr) - /2^g (^s, y) + {Eo + M' + Ag+) cot(s - ^o) 

ni n2 

- (Eo + Mg+ - M- + Ag-) cot(s + ^0) - E ^ij"^^ cot(s - s^f )) - £ A^if ^ tan(s - sf^^). (6.33) 

To recover the reflected, surface and diffracted waves, we also need to continue the function 
Se{s) to the strip — 37r/2 < Res < 0. In the strip — vr < Res < 0, the function Fie(s) has poles 
at the points — ^O) — '^ + ^0 (the geometric optics poles) and may have surface wave poles at the 
zeros of d{kQCOss) in the domain — vr < Res < 0, Im(A:o cos s) < and at the zeros of d{kQSva.s) in 
the strip — vr < Res < 0. Denote the residues of Fie(s) at the surface wave poles, s);- , as A'^g- 
(j = 1, . . . , m'^). This brings us to 

S)^'{s) = -— ^ / Fie(o- - vr) cot((T - s)da + A^ cot(s - 6*0) 

ZTII J —ioci 

- A+ cot(s + 0o) + E A^i]^^ cot(s - s^J.-^^), -vr < Re s < 0. (6.34) 

i=i 

By shifting the contour of integration further to the left by vr/2 we find the analytical continuation 
of the function Se (s) into the strip — 3vr/2 < Res < vr. In doing so we take into account the poles 
of -Fie(s) at s = — vr — 9q and Sg • with the residues A'^g ■ [j = m[ + 1, . . . , mi). 

The analytical continuation of Se2is) from the strip —vr/2 < Res < vr/2 into the strip — 3vr/2 < 
Res < —vr/2 is implemented by shifting the contour of integration in ()6.13p to the left by vr. The 
poles of F2e{s) in the strip — vr < Res < comprise the the geometric optics poles —vr/2 + 9q and 

(2L\ (2L) 

—vr/2 — 9q and the surface wave poles Sg- with the residues A'^g- (j = 1, . . . ,1112) at the zeros of 
d{kQCoss) in the domain — vr < Res < 0, Im(A;o cos s) < and at the zeros of d{kQsms) in the strip 
— vr < Res < 0. This enables us to write the analytical continuation of the function Se{s) into the 
strip — 3vr/2 < Res < as 

5e(s) = /L(s, -^) + iLis, 0) + (£^0 + A-) cot(s - ^o) 

-(^0 + A^) cot(s + ^o) + E ^if ^ cot(s - sif ^), -^ < Res < 0, 
5e(s) = Ii=g(s, -vr) + /*g(s, -tt) + (^0 + A^ + M+) cot(s - ^o) - (^0 + A+ + Mg") cot(s + ^o) 



vr 



+ E <^ -Otis - sl-^) + E Kr *--(« - <T^ -vr < Re s < -|, 
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Se{s) = -II ( s, -^) + Ilis, -vr) + {Eo + A- + M+) cot(s - ^o) 
-(So + A+ - A" + M-) cot(s + 0o) 



+ Y: N^r ^°t(« - 4i ^) + E <^ tan(s - sif ^), -:^ < Re . < -vr. (6.35) 

Here, s^- (j = 1, . . . ,m'() and Sg- (j = 1, . . . ,'"1.2) are the surface wave poles of the functions 
Fie(s) and F2e{s), respectively, in the domain —tt/2 < Res < 0. 

In the case of the function 5/j(s) all the formulas derived can still be used if we replace /3 by vr — /3 
and interchange the parameters 7^ with 74 , the functions <1>^ (s), "^i with ^2 (-s)) ^2 ' respectively. 
The constants A^, and M^ corresponding to the function Sh{s) are denoted by A^ and M^ . 

6.2 Incident, reflected, surface and diffracted waves 

Let C be a closed contour that comprises the Sommerfeld contours 7+ and T- and two contours 
Ql and Gr. The starting points of the contours Ql and Gr are —3-k/2 + ioo and 37r/2 — ioo, and 
they pass through the points s = — vr and s = vr, respectively. The contour Gji is described by the 
equation Res = vr — gd(Ims) sgn(Ims), and the contour Gl is symmetric to Gr with respect to the 
origin. Here, gdx is the Gudermann function gdx = cos~^secha;. As k^r — )• 00, the electric and 
magnetic field can be represented as 

E^^El+El + Et + Ei 

H^^ Hi + HI + HI + Hi, kr -^ 00, (6.36) 

where El, HI are the incident waves, E^, H^ are the reflected waves, Ef,, H^ are the surface waves, 
and E'i, H^ are the diffracted waves. 

Write down first those poles of the functions Se{s) and Sh{s) which are in the interval (— vr,vr) 
and which generate the incident and reflected waves. These poles are 

s = -e + eo£{-TT/2,Tr/2), s + 0e (0,vr/2), 

s = vr-6'-6lo e (0,7r), s + 6* G (7r/2, vr), 

s = -e-eo£{-TT,0), s + 0e (-vr/2,0), 

s = Tr-e + eoe {tt/2, vr) if 6^0 < 61 < vr/2, s + 6* G (vr, 3vr/2), 

s = -vr - 6* + 6*0 G (-7r, -vr/2) if < 6* < 6^0, s + 6* G (-vr, -vr/2). (6.37) 

The incident and reflected electrical and magnetic waves can be recovered from (j6.ip . ()6.3ip . (j6.32p 
and (J635]) 

j^i _ ^ -ifcopcos{e-6»o) jji _ ^ g-ifcopcos(e-eo) 

El = -{Eo + A+)e-^'=o^=°'^(^+^o) - {Eo + M+)e^^'«''™" (''+''«) 
+uj{e; 0, 9o){Eo + A- + M+)e^'^°^™^(^-^°) + w(0; Oq, vr/2)(So + A+ + M-)e'^'^P^°<^-^°\ 
ZHl = -{ZHo + A+)e-*'=«^"°'^(^+^o) - {ZHq + M+)e*'^«''"°^(^+^o) 
+ w(6';O,0o)(^^o + A;^ + M+)e*'=«^"°'^(^~^o)+w(0;0o,vr/2)(^^o + A+ + M^-)e^^o^™"(^~''"), (6.38) 
where 
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To recover the surface waves we select those surface wave poles which meet the conditions 



< Re(- 



ej > 



<7r-5 ', J = 1,2, 



.ni<ni, 1 = 1,2, 



^-5]'''^<Re(-0 + 4''^)<O, j = l,2,...,mi<mi, 1 = 1,2, 



.40) 



where g, = gd(ImSg- )sgn(ImSg- ), P = R,L. Denote these surface wave poles as Sg- 
(ReSg- > 0,j = l,...,hi) and Sg- (ReSg- < 0,j = l,...,m/), and the associated residues 
as TVg- and N^- , respectively. Let 



MR) 



(IB) 



,(IR) 



m 



ML) 



~SlL)^ 



X):p' = max{0, -TT + gp + Re(st/0}, Kj = min{7r/2, tt + g^ + Re{s)!j>)}. 



(6.41) 



We note that the poles Sg • and Sg • generate the surface waves if and only if —-K+g- +Re(Sg ■ ) < 

7r/2 and vr + gf^'^ + Re(sJ:^^) > 0, respectively. By making use of ^^, (fOB . (f02D and (lOSjl we 
write down the surface waves Et and Ht as 



K = Y. 



1=1 



ni 



- Y: KT^{e, Aif , 7r/2) exp{-ifcopcos(0 - sif )} 



m; 



+ Y. K;^u;{e, 0, Ai;^^ exp{-ikopcos{e - s)!j'^)} 






"i 



,(ZR) 



Y N),Puo{e, A},'f , 7r/2) exp{-ifcopcos(0 






)} 



m; 



i=i 



,(«/?) 






(6.42) 



(IR) 



where hi, rhi, AL- and AL- are the counterparts of the corresponding numbers hi, hii, X- and 

Ag • in the magnetic case. 

The diffracted waves are determined in a standard manner in the form 



E' 



pikop pikop 



/k)P 



De{9), ZHl 



'h)P 



Dh{0), 



(6.43) 



where De{9) and Dh{9) are the diffraction coefficients 



De{e) 



-i-k/A 



\Se 



7r)-5e(e + 7r)], Dh{e) 



-i-w/A 



\Sh{e-^)-Sh{e + ^)]. 



(6.44) 



Conclusions 

We have found a closed-form solution of the classical problem on scattering of an electromagnetic 
plane obliquely incident wave from an impedance wedge in the right-angled concave case. By split- 
ting the incident wave into two waves, the one that strikes the vertical wall and the other that 
impinges upon the horizontal face of the wedge, we have reduced the boundary value problem for 
the Helmholtz equation to two vector RHPs subject to the symmetry condition for the unknown 
vectors, $+(77) = <1>~(— r/). The matrix factorization problem has been solved by the technique of 
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the RHP on a genus-3 Riemann surface in the case of the boundary conditions p.l|) . In the general 
case (jl.2p of the boundary conditions the matrix can be factorized as well. However, this requires 
tedious computations of the analogues of formulas (I3.22|) . ()4.ip and ()4.2p . Due to the symmetry 
of the problem we have managed to avoid the genus-3 Riemann ^-function and to solve the Jacobi 
inversion problem in terms of elliptic functions. The Wiener-Hopf matrix- factors have been found 
by quadratures and, eventually, a closed-form solution to the two vector RHPs associated with the 
diffraction problem of interest has been constructed. In addition, we have proved the uniqueness 
theorem (Theorem 5.1) which generalizes the known result in the scalar case (see for example (J18j) ) 
to the vectorial case. The solution is unique when all the four zeros of the characteristic polynomial 
(I5.65P of the problem lie in the upper half-plane and has n arbitrary constants with n being the 
number of the zeros lying in the lower half-plane. In the uniqueness case we have shown how to 
express the spectra of the problem and therefore the Sommerfeld integrals through the solution of 
the two vector RHPs. By using these expressions we have analytically continued the spectra to the 
left and to the right and by steepest descent method recovered the incident, reflected, surface and 
diffracted waves. 
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